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ABSTRACT 

Context. Some rapidly variable extra-galactic radio sources show very high brightness temperatures Tg > 10 12 K and high degrees of circular 
polarisation (~ 1%). Standard synchrotron models that assume a power-law electron distribution cannot produce such high temperatures and 
have much lower degrees of intrinsic circular polarisation. 

Aims. We examine the synchrotron and inverse Compton radiation from a monoenergetic electron distribution and discuss the constraints placed 
upon it by radio, optical and hard X-ray/gamma-ray observations. 

Methods. The standard expressions of synchrotron theory are used. Observational constraints on the source parameters are found by formulating 
the results as functions of the source size, Doppler boosting factor, optical depth to synchrotron self-absorption and maximum frequency of 
synchrotron emission, together with a parameter governing the strength of the inverse Compton radiation. 

Results. The model gives brightness temperatures Tg ~ 10 13 to 10 14 K for moderate (< 10) Doppler boosting factors together with intrinsic 
degrees of circular polarisation at the percent level. It predicts a spectrum /,, oc y 1/3 between the radio and the infra-red as well as emission in 
the MeV to GeV range. If the energy density in relativistic particles is comparable to or greater than the magnetic energy density, we show that 
electrons do not cool within the source, enabling the GHz emission to emerge without absorption and the potentially catastrophic energy losses 
by inverse Compton scattering to be avoided. Magnetically dominated sources can also fulfil these requirements at the cost of a slightly lower 
limit on the brightness temperature. 

Conclusions. We suggest that sources such as PKS 1519 -273, PKS 0405 -385 and J 1819 +3845 can be understood within this scenario 
without invoking high Doppler boosting factors, coherent emission mechanisms, or the dominance of proton synchrotron radiation. 
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1. Introduction 

The rapidly varying radio flux density observed in several 
extra-galactic sources implie s a very high brigh tness tempera¬ 
ture 7 b at the source JWagner & Witze]lfl9 95). In cases such 
as PKS 1519 -273 and PKS 0405 -385, where the vari¬ 
ability is identified as inte rstellar scintillation iMacauart et alJ 
2000; iRickett et al .1120021) realistic models of the scattering 
screen require 7 ’b > 10 13 K. The recently discovered diffrac¬ 
tive scintillation in J 1819 +3845 requires T B > 2 x 10 14 K 
dMacauart & de Bruvnll2005i) . In other sources, the variability 
can be interpreted as intrinsic, in which case the implied tem¬ 
perature can be even higher. 

Such high temperatures are difficult to understand within 
standard models of synchrotron emission, which assume 
a power-law distribution of electrons. The energy radi¬ 
ated by inverse Compton scattering in the Thomson regime 
rises dramatically when the intrinsic brightness tempera¬ 


ture exceeds a certain threshold, roughly equal to 3 x 
10 n K in the case of a power-law electron distribution, al¬ 
thoug h a flaring source might exceed this limit for a s hort 


time 

iKellermann & Paulinv-Toth 1969:|sivsh 1992): Melrose 

bnni 

. However, as first pointed out by Crusius-Waetze 


( 1991 ). the threshold temperature is higher for a monoenergetic 
electron distribution than for a power-law distribution, since, 
in the former case, photons of low frequency (compared to the 
characteristic synchrotron frequency) can escape without ab¬ 
sorption by low energy electrons. 

High brightness temperature sources frequently display 
circular polarisation at the 1% level or above dMacauartl 
booi) . Because this is much higher than the value m e c 2 /k B T B 
conventionally estimated for the intrinsic emission of a 


power-law electrc 
favoured explana 

>n distribution, propat 

ration effects are the 

tion (Jones & O’Dell 

1977b; Wardle et al. 

1998; Macquart & Melrose 2000; Beckert & Falcke 

2002; 

Ruszkowski & Be 

gelman 2002; Broderick & Blandford 

2003; 

Wardle & Homan 

2003). A monoenergetic electron model, on 


the other hand, predicts intrinsic circular polarisation at the % 
level, obviating the need for a conversion process. Additional 
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predictions are a hard I Y oc y 1/3 synchrotron spectrum extend¬ 
ing from the radio at least up to the infra red, and an inverse 
Compton component in the MeV to GeV range. 


2. Model parameters 

We consider an idealised model characterised by a single length 
scale R and Doppler boosting factor D = y/l -/3 2 /( 1 —/J cos </>), 
where c(3 is the source speed with respect to the rest frame of 
the host galaxy and cp is the angle between the velocity and the 
line of sight. A homogeneous, isotropic, monoenergetic dis¬ 
tribution of electrons with Lorentz factor y and number den¬ 
sity n, embedded in a uniform magnetic field B is assumed in 
the rest frame of the source. This captures the relevant prop¬ 
erties of a power-law electron distribution with a lower en¬ 
ergy “cut-off”, provided the synchrotron opacity at low fre¬ 
quency is dominated by electrons whose energy is close to 
that of the cut-off, as is the case if the distribution rises to¬ 
wards the cut-off sufficiently rapidly: dlnn/dlny > -1/3. 
Quasi-monoenergetic distributions of this type can account for 
the lack of Faraday d epolarisat ion in parsec- scale emission 
regions I' Wardlei 119771 Jones & O’Delll Il977al) and have re¬ 
cently been discussed in connection with statistical trends in 
the observed distribution of superluminal velocitie s as a func- 
tionof observing frequency and redshift (Gon al-Krishna et all 
2004). Monoenergetic models have been considered previ- 
ousl y in connection with h igh brightness temperatu re sources 
bv jCrusius-Waetzeil dl99lh . and, more recently, by Prot heroel 
(2003), who, however, adopted the assumption of equipartition 
between particle and magnetic energy densities in the source 
and introduced an additional, potentially large, geometrical fac- 
tor. I.SIvsh! (1 992! ) also treated a monoenergetic model, but did 
not allow for the possibility of multiple Compton scatterings 
and further restricted his treatment to optically thick sources. 

Neglecting, for a moment, the angle between the line of 
sight and the magnetic field, the source model contains five free 
parameters: R, D, y, B and n. Of these, only two are directly 
constrained by observations. Firstly, an upper limit on the linear 
dimension R, is given by the maximum permitted angular size 
inferred from the presence of scintillation. Secondly, surveys of 
superluminal motion ( iCohen et al J2003t) reveal apparent trans¬ 
verse velocities that are mostly less than 10c, but extend up to 
30c. These data suggest that D < 10 for most sources, although 
it may be up to 30 or 40 in a few cases. If the source variability 
arises intrinsically, rather than from scintillation, an upper limit 
to the quantity R/D 2 is obtained. 

Additional constraints can be found by eliminating y. B and 
n in favour of three new parameters. The first of these is the 
optical depth to synchrotron self-absorption r s . Denoting the 
Thomson optical depth of the monoenergetic electrons by tt = 
ncrjR (where cr T is the Thomson cross-section) one has 


T s = 


V3 r T m e c 3 K 5 /i (x) 


( 1 ) 


87re 2 v s y 3 

(e.g., tLongaidl992l) ). Here, Ksp(x) is a modified Bessel func¬ 


tion. 


( 2 ) 


v is the observing frequency, v s = 3eB sin 6y 2 /(4nm e c) is the 
characteristic synchrotron frequency, for an angle 9 between 
the line of sight and the magnetic field B (both measured in the 
frame in which the source is at rest) and z is the redshift of the 
host galaxy. t s is a convenient parameter because it controls the 
brightness temperature, which has a single maximum close to 
r s = 1, when the other parameters are fixed. 

The second new parameter, £, determines the inverse 
Compton luminosity. Synchrotron photons act as targets for in¬ 
verse scattering by the same electron population that produced 
them, giving rise to a first generation of upscattered photons. 
These, in turn, act as targets off which the electrons produce 
a second generation, and so forth. We define as the ratio of 
the energy densities (or, equivalently, luminosities) in consecu¬ 
tive generations, assuming the scattering processes occur in the 
Thomson regime (e.g., lMelrosell2002l) . This quantity is some¬ 
what sensitive to the geometry and homogeneity of the source 
(jProtheroel2()02i) . but, for a uniform, roughly spherical source 
that is optically thin to scattering, one can write 

? = 4y 2 r T /3 (3) 

since the photon frequency is increased on average by the fac¬ 
tor 4y 2 /3 per scattering, and the probability of scattering is tt. 
The assumption of Thomson scattering leads to a divergent lu¬ 
minosity whenever f > 1, independent of the source of the ini¬ 
tial (zeroth generation) photon targets. This phenomenon has 
acquired the name Compton Catastrophe. If the source is opti¬ 
cally thin to most of the synchrotron photons, £■' is equal to the 
ratio of the energy density in these photons to that in the mag¬ 
netic field. This is, in fact, the most commonly used definition. 
However, using our definition 0 the condition for catastro¬ 
phe remains £ > 1 in both the optically thin and thick (to syn¬ 
chrotron absorption) cases as well as in the case where photons 
of the cosmic microwave background provide a more effective 
target than do the synchrotron photons. We show below that in¬ 
verse Compton scattering indeed takes place at least initially in 
the Thomson regime. Therefore, in order to avoid inordinately 
large Compton losses, we require 

£ < 1 (4) 


which automatically ensures tj <s 1. 

The observed brightness temperature 7 b, related to the spe¬ 
cific intensity of radiation /,, in the direction of the source by 
7 b = c 2 l v /( 2 v 2 k H ). follows straightforwardly from the solution 
of the equation of radiation transport: 


^Tb = D I yF(x) \ 
m e c 2 1 + z \2x 2 K 5/3 (x)) 


(5) 


poo 

where F(x) = x J dtKs 13 (f) is the standard synchrotron func¬ 
tion in the Airy integral approximation. Eliminating the param¬ 
eters tt, v s and y using Eqs. 0. 0 and 0 . gives 


m e c 2 


/3 3/ V\ 1/5 / CD 6 \' /5 /l - e~ Ts 

\ 4 5 7re 2 v j \(l+z) 6 j \ r s 1/5 


F(x) 

, x9I5K 5/3 W 


( 6 ) 


The first term in parentheses on the right-hand side of this 
equation is independent of the source parameters. The second 


x - v(l +z)/(Dv s ) 
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is constrained by observation to be < 10. The third reaches a 
maximum of the order of unity at r s ~ 1. The fourth, how¬ 
ever, diverges for small x as x~ 2/15 . Thus, even with £ < 1 
and D < 10 , it is possible to choose an x for which this for¬ 
mula gives an arbitrarily high brightness temperature at any 
specified observing frequency. However, this model predicts 
a synchrotron spectrum that rises with frequency at least as 
fast as I v oc v 1 ' 3 between the observing frequency v and the 
frequency v/x. In principle, this can be constrained by obser¬ 
vation. For PKS 1519 -273 and PKS 0405 -385, for exam¬ 


ple, the modest o ptical fluxes suggest that v max < 10 14 Hz 
IHeidt & Wagnei 119961 and Wagner, priv. comm.), although 


it is not known whether these observations coincided with 
an episode of high brightness temperature radio emission in 
these variable sources. Nevertheless, for our third parameter, 
we adopt the frequency v max = v/x above which the optically 
thin synchrotron radiation cuts off. 


3. Brightness temperature and circular 
polarisation 

Rewriting the brightness temperature Eq. <01 in terms of the 
three new parameters r s , £ and v max , expressed in convenient 
units, we find that it is independent of the source size, and de¬ 
pends only weakly on £ and v max : 


T b = 1.2 x 10 


14 


£) 6 


1/5 


Kri 


(1+z) 6 


1 - e“ 

1/5 


2/15 —1/3 tv- 

V maxl4 V GHz 


(7) 


where v ma x !4 = v ma x/(10 14 Hz), v GHz = v/(l GHz), D ,0 = 
D/ 10, and the approximations Fix) ~ 2.15X 1 ' 3 and K^p(x) ~ 
1,43x 533 , valid forx «: 1, are used. According to this equation, 
brightness temperatures of up to roughly 10 13 K at GHz fre¬ 
quencies, as observed in PKS 1519 -273 and PKS 0405 -385, 
can be achieved with £ < 1 and 'D - 1. With D ~ 15, the 
model permits Tb = 2x10 14 K, as observed in J 1819 +3845. 

In addition to the high brightness temperature, a particu¬ 
larly interesting source property is the degree of intrinsic cir¬ 
cular polarisation r c . Assuming a pure electron-proton plasma 
that is optically thin to synchrotron self-absorption, this quan¬ 
tity is also independent of source size: 


—f- 

3y \x 
= 0.019 x 


1/3 


cot6»r(l/3) 


(1 +z)t s 
£>io£ 


1/5 


v 1/5 . . COtP 
max 14 


( 8 ) 


iMelrosel 1980 ). In the case of a power-law electron distribu¬ 
tion, r c changes sign when the optically thick regime is entered 
Ijones & Q’Delill977bh . We are not aware of the correspond¬ 
ing calculation for a monoenergetic distriubtion, but, to order of 
magnitude, one can estimate the peak value using this expres¬ 
sion, which is remarkably insensitive to all source parameters 
other than the magnetic field direction. Several extra-galactic 
sources of extremely high brightness temperature display cir¬ 
cular polarisation at the percent level iMacauartl2003l) . in par¬ 
ticular PKS 1519 -273 and PKS 0405 -385. In the absence 
of a low-energy cut off in the electron distribution, r c ~ l/y. 


which is far too small to explain the observations. However, 
Eq. <0 shows that for a mono-energetic electron distribution, 
the intrinsic emission can be polarised at the % level, or above, 
depending on the geometry of the magnetic field configuration. 


4. Discussion 


The electron Lorentz factor implied by the above analysis is: 


y = 2.8 x 10 3 




,2/3 V 

max 14 


+s(l +Z) 


1/5 


w 1/3 

y GHz 


(9) 


A key ingredient of this model is the absence of electrons 
of lower Lorentz factor, since these would absorb the GHz 
emission, leading to a reduction of the brightness tempera¬ 
ture. Specifically, we require a quasi-monoenergetic distribu¬ 
tion such that d In n/d In y > -1/3 at Lorentz factors lower 
than that given by Eq. < 0 . Such a distribution is not a natu¬ 
ral consequence of, for ex amp le, the first-order Fermi process 
at relativistic shocks (e.c.. lKirkl2005l) . On the other hand, a rel¬ 
ativistic thermal distribution, which rises at low energy as y 2 is 
well-approximated by a monoenergetic distribution of energy 
roughly equal to the temperature. The addition of a power-law 
tail to higher energy would not change this conclusion. 

The Lorentz factor implied by Eq. 0 is higher than the 
cut-o ff conventionally assumed when modelling radio sources 
(e.p. JGopal-Krishna et al.fl2004i) . Nevertheless, scenarios exist 
which suggest such values. One example is an electron-proton 
jet with a bulk Lorentz factor T ~ 10 which is thermalised at a 
shock front. If the downstream electron and ion temperatures 
are equal, the distribution can be approximated as monoen¬ 
ergetic with an electron Lorentz factor of T m p /m e , where m v 
and m e are the proton and electron masses. Another possibil¬ 
ity is that the electrons are accelerated at a current sheet in an 
electron-proton plasma in which the magnetic energy density is 
comparable to the rest-mass energy density dKirl 12004 . Each 
of these possibilities relies on the composition of the source 
plasma being electron-proton. Interestingly, so does the rela¬ 
tively high degree of intrinsic circular polarisation given by 
Eq. 0. 

Although it is conceivable that continuous re-acceleration 
prevents the accumulation of low energy electrons, both the 
current sheet and the shock scenario envisage a finite escape 
rate of particles from the acceleration or thermalisation region. 
Escaping particles subsequently cool by synchrotron and in¬ 
verse Compton emission. Therefore, the model electron distri¬ 
bution is self-consistent only if these particles can be evacuated 
from the source in a time shorter than the cooling timescale. 
The ratio of the electron cooling time / c „ 0 | to the light-crossing 
time of the source can be written as: 


ctcooi/R = nit; 


GO) 


where 77 is the ratio of the energy density in relativistic electrons 
to that in the magnetic field. Writing R = 0.01 R -2 parsec, we 
find: 


V = 


ynmc 


(B 2 / 8 tt) 


= 2.9 


D! 3 * 8 


10 ^ 


(1 +z) 13 r 3 v 8 I4 


1/5 


v G L^2 sin 2 0(11) 
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Clearly, very small sources tend to be particle dominated and, 
since £ < 1, they satisfy the self-consistency requirement 
ctcool/R > 1- However, Eq. CD shows that r\ is also quite sen¬ 
sitive to £>, and v max , so that magnetically dominated sources 
are by no means ruled out, provided they have ^ < r\. Since 
the brightness temperature is proportional to (Eq. 0 it is 
slightly lower for magnetically dominated sources. 

Although lacking a compelling physical justification, the 
assumption of equipartition, r\ = 1, can be used to define 
an “equipartition Doppler factor”. This leads, in the standard 
model, to a relatively low limit on the brightness temperature 
Tr Pn < 3 X 10 10 K Isingal & Gopal-KrishnaifToSSi: iReadheadl 
1994 ), which has some observational support ( Cohen e t al. 
2003). In the monoenergetic model, however, this assumption 
re-introduces a dependence on the source size, but does not 
substantially constrain the brightness temperature, as can be 
seen from equations 0 and ill It . 

The parameter 


yhv s 
m e c 2 


2.3 x 10^ 4 


(1 + z) 4 ^ 


v 17/t5 -1/3 
max 14 uHz 


( 12 ) 


which gives the ratio of the energy of a photon of the charac¬ 
teristic synchrotron frequency to the electron rest-mass, as seen 
in the rest frame of a relativistic electron is also independent of 
source size. For £ <s 1, the first inverse Compton scattering 
takes place in the Thomson regime. In this case, it is consis¬ 
tent to require £ < 1 in order to avoid an excessively large en¬ 
ergy demand on the source i.e, in order to avoid the Compton 
Catastrophe. The first generation of inverse Compton photons 
has a frequency of approximately 


v l 


4.3 


D 2 (P- 


(1 + z) 2 t 2 


1/5 


v 19/1 , 5 4 VrH /3 MeV 


(13) 


and its flux can be estimated to be 


Tic ~ 4.5 x 10“ 


(1 +z)W 


TP- 


10 


1/5 


1/15 1/3 fjr 

V maxl4 V GH Z 7"GHz 


(14) 


where 'Fqw, is the flux observed in the radio at frequency 
vghz GHz. This estimate of the inverse Compton flux is 
generally above the detection threshold of instrument^ on 
the INTEGRAL satellite, as noted by iProtheroel (12003}) . 
Subsequent generations of inverse Compton scattered photons 
are likely to fall into the Klein-Nishina regime, and the max¬ 
imum photon energy achieved by multiple inverse Compton 
scattering is ultimately limited by the electron energy, as seen 
in the observer’s frame, which takes the value: 


Dymc 2 

1 + z 


14 


^io£ 

(1 +z) 6 t s 


\ 1/5 


y 2/15 -1/3 

max 14 uHz 


GeV 


(15) 


To summarise, synchrotron radiation from a monoenergetic 
electron distribution reproduces the extremely high bright¬ 
ness temperatures observed in variable extra-galactic radio 
sources, and explains the observed levels of circular polarisa¬ 
tion. Therefore, it d oes not appear necessary to appeal to co¬ 
herent mechanism s llKrish nan & Wiitall990tlBenford & Leschl 
1 998 : Begelman et a 1.120051) or to proton synchrotron radiation 


iKardashevfooQOt) to understand these objects. Testable predic¬ 
tions of the theory are a hard radio to infra-red spectrum and 
gamma-ray emission in the MeV to GeV range. 
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